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Orbital Angular Momentum (OAM) [1] of light has 
recently attracted a lot of attention as a new promis¬ 
ing resource for fundamental and applied physics [2-10], 
such as biophysics [2], microscopy [3], astronomy [4] and 
metrology [7]. In this framework, Circular-Beams (CiBs) 
were introduced in [11] as a very general solution car¬ 
rying OAM of the paraxial wave equation. The impor¬ 
tance of the CiBs is underlined by the fact that they 
represent a generalization of many well known beams 
with OAM. Indeed, for particular values of the beam 
parameters, the CiBs reduce to standard [12] or elegant 
Laguerre-Gauss [13], Hypergeometric-Gauss [14, 15], Hy¬ 
pergeometric [16] modes and others [11]. Then, the study 
of the CiB properties allows to better understand, in a 
unified way, the features of such well known beams. For 
instance, since the propagation of CiBs through ABCD 
optical systems can be easily described in terms of their 
defining parameters qo and q± [11], the same propagation 
law can be also applied for their known particular cases. 
Moreover, while Laguerre-Gauss (LG) beams seem a nat¬ 
ural choice for describing OAM states in relatively simple 
terms, in several cases they only represent an approxima¬ 
tion of the generated beam. Indeed, a special CiB quite 
different from an LG mode, is obtained by applying a sin¬ 
gular phase factor exp (imcft) to a Gaussian beam (see eq. 
(18) and [14]). This technique represents a very common 
and simple way to generate OAM light [5, 8, 10]. Thus, 
the study of the Circular-Beams gives insight for the pre¬ 
cise modeling, generation and exploitation of OAM in 
classical or quantum applications. 

In this letter, we analyze the properties of the CiBs: 
in particular, we will explicitly derive their expansion 
in terms of Laguerre-Gauss modes and their normaliza¬ 
tion constant, allowing the evaluation of their free-space 
divergence. We will also investigate the constraints on 
the values of the beam parameters arising by the request 
of square integrability and regularity at any finite point. 
Then, we will study the experimental generation of a sub¬ 
class of CiBs. Our results will provide an useful tool for 
the complete modeling of CiBs, paving the way towards 
their experimental implementation. 

CiBs were introduced in [11], without the evaluation 
of their normalization. As a first result of this letter, we 
give the detailed expression of the normalized CiB. Let’s 
define the complex parameters q{z), q{z), £ and the scale 


factor x(z) as: 


q{z) = z + q 0 , 


q{z) = z + qi. 
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were the real and imaginary part of the two complex 
beam parameter qo and q± are given by qo = izo — do, 
qi = iz\ — di being zq = ttWq /A and do the analogues of 
the confocal parameter and the location of the waist for 
a Gaussian beam. The normalized monochromatic CiB 
in cylindrical coordinates x = (r, <j>, z) and propagating 
along the z axis is defined as exp[i(wt — /cz)]CiB^ 91 ^ (x) 
where 
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In the above equation \F\ is the confluent Hypergeomet¬ 
ric function iFi(a, 6; a) = Et=o Kyr'ct+nj AT and "4.™ 
a normalization factor given by 


*(« r(n-e)r(n-£)|m|! 2n 

h r (-I) r (-^)n!(H+n)! l?l 

= 2 Fi[—|,- 4 l + |m|,|£| 2 ]. (3) 

The parameters m £ Z and p € C respectively rep¬ 
resent the amount of OAM carried by the beam and 
the analogous of the Laguerre-Gauss mode radial index. 
As noted in [11], the field is invariant (up to a nor¬ 
malization factor) by the transformation (p,m,qo,qi) —► 
(—p — 2|m| — 2,m,qi,qo) 1 - As we will see, to guarantee 
the square integrability and regularity at each finite z 


1 In [11] the field was defined with p replaced by i'y — |m| — 1 
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and r, some constraints will arise on the allowed values 
of p. 

As a second result of this letter, we give the expan¬ 
sion of the CiB in terms of the Laguerre-Gauss beams 
LG njm (x) when q\ ^ q$ and 9m(go) > 0: 


+ OO 

CiB^(x) = Y, A^LG^r, </>, z - d 0 ), (4) 

n =0 

where 
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The above expansion does not hold when q\ = g(j 
and 9m(go) > 0: indeed, in this case, the LG„ im 
mode is directly obtained as a particular CiB, namely 
LG H:m (r, (/), z - d 0 ) = (-l) n CiB^°£ o) (x). Equation (4) 
implicitly shows that the CiBs satisfy the paraxial wave 
equation since they can be represented as a linear com¬ 
bination of its solutions (i.e. the LG modes). For com¬ 
pleteness, we here report the explicit expression of the 
standard Laguerre-Gauss modes: 
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where n, to G Z with n > 0, r) is the generalized 

Laguerre polynomial, W(z) = W 0 yjl + (z/z 0 ) 2 the beam 
size and exp[i£(z)] = (zq + *z)/|zo +iz\ the Gouy phase. 

We now demonstrate eqs. (3) and (4). When qi ^ q$ 
and 9m(go) > 0 we may exploit the expansion of the Hy¬ 
pergeometric function in terms of Laguerre polynomials: 


HZ y? r(o + n)T(b) Y n L { n~ 1] ( Z) 
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Indeed, by using Z = w ^ do) 2 and Y = 
in the above equation, relation (4) can be easily ob¬ 
tained. Regarding the value of the constant (see 

eq. (3)), it can be determined directly from the expansion 
given in eq. (4). By recalling that the Laguerre-Gauss 
modes are orthonormal with respect to the scalar prod- 
uct J dcj) J dr r LG nirrai (x)LG„ 2jm2 (x) d ni ^ n2 S rni ^ rn21 
the CiB is easily found to be normalized to 1 if the con¬ 
stant is defined as (3). 

As already anticipated, to obtain square integrable 
and regular beams, some constraints (summarized in eq. 
(8)) arises on the allowed values of p. Let’s start by 
the square integrability, equivalent to the convergence 
of the sum in eq. (3). The latter converges to the 
Hypergeometric function 2 Ei in three cases: a) when 
9m(g 0 ) > 0, 9m(gi) > 0, Vp; b) when 9m(g 0 ) > 0, 


9m(gi) = 0,+oo and 5ie(p) > —1 — \m\; c) when 
9m(go) > 0, 9m(gi) < 0 and p = 21 with £ G N. In 
the latter case the sum becomes finite, since in eq. (4) 
the expression 5^+f^r(n — |)/T(—|) must be replaced 
by n)!. The three cases above cited 

respectively correspond to |£| < 1, |£| = 1 and |£| > 1. 
As noticed in [11], the same constraints can be obtained 
by looking at the CiB behavior at large r. 

However, further conditions on p arise by requiring the 
regularity of the field at each finite z. In the cases a) 
and c) the field is regular at any finite r and z. On 
the other hand, when 9m(go) > 0 and 9m(gi) = 0 the 
field may be singular at z = d±. Indeed, since in this 
case we have lim z _ >dl CiB^ l?1 ^(x) oc r p+ l m l, regularity 
in r = 0 requires 5fte(p) > — \m\. Then, the integrability 
and the absence of singularity at any finite r and z are 
only guaranteed in the following four inequivalent cases: 

I) 3m(g 0 ) > 0, 9m(gi) = +oo , Sfte(p) > — |m| — 1 

II) 3m(g 0 ) > 0, 9m(gi) < 0, p = 2£ with t G N 

III) 3m(g 0 ) > 0, 9m(gi) = 0, Ute(p) > — \m\ 

IV) 3m(g 0 ) > 0 , 9m(gi) > 0 , Vp G C (8) 


Other equivalent cases may be obtained by the symmetry 
(Pi m, g 0 , gi) (-p- 2|m| - 2,m,gi,g 0 ). 

To further investigate the properties of the CiBs, we 
now evaluate their free-space divergence, an important 
parameter for experimental implementation and propa¬ 
gation through long distances. For general light beam, 
the root mean square (rms) divergence 0 rms can be de¬ 
fined as [17-19] 

0 rms = lim (9) 

Z—¥+ OO Z 


where cr^ ms (z) is the second moment of the intensity /(x): 

p2tt p+oo 

iW = / # / drr [r 2 /(x)] . (10) 
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By exploting the expansion of the CiB in 
terms of the LG modes (4) and the relation 
2 /d<^ fdrr 3 LGq m (x)LG* m (x) = W 2 (z)[B m>e 6 n , t - 
C m ,eS n + i,e - Cm in S n ,i+ 1 ] with = |m| + 2n + 1 and 

Cm,n{z) = sjn{\m\ + n) exp[2i£(z)] it is possible to give 
an explicit expression for of ms (z), namely: 

2a? W 2 (z) d0) = 1 + H + + Ke[X( " )] ’ (U) 

where X{z) = £(p — 4>p%) exp[2i£(z)] and 
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By plugging the above result into (9), the rms divergence 
of the CiBs may be expressed as 
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FIG. 1. Divergences 9 lms and $ee in function of \m\ for the 
CiB with qi = Ke(qo) and p = 0, p = — |m| + 1 or p = — \m\. 
When p = —|m|, only $ee can be defined. The beam is only 
defined for integer values of m and lines represent just a guide 
for the overall trends. Note that when p = 0 the CiB p , m 
reduces to the Laguerre-Gauss mode LGo, m - 


In the above equation 9 0 = Wq/(\/2 zq) is the divergence 
of a Gaussian Beam (i.e. a CiB with p = m = 0) with 
confocal parameter zq and beam waist Wq = \JXzq/tt. 
As we will show in the following (see eq. (17)), the 
divergence 9 lms is not finite when = 0 and 

3?e(p) = — |m|. Indeed, for such beams, at fixed 0 and 
large r, the intensity fall-off as r -4 . While the field is 
thus square integrable, the integral (10) is divergent and 
the <Trms(z) cannot be defined. 

To avoid such problem, we may use an alternative 
definition of the divergence in terms of the so-called 
encircled-energy. It requires the determination of the 
radius Ree(z) whose corresponding circle centered on 
the beam axis contains a given amount of beam energy. 
Then, -Ree(-z) must be calculated by the implicit rela¬ 
tion f d (f> J^ EE ^ d rrl(x) = I 0 , with I 0 a fixed constant. 
Given Ree(z), the corresponding divergence $ee can be 
determined similarly to eq. (13): 

»EE= lim 5sA>. (14) 

z— H-oo Z 


If we choose Io = 1 — 1/e (chosen to achieve 9ee = 9 rms = 
9 0 for the Gaussian beam), the divergence 9ee of a CiB 
can be obtained by the following implicit relation: 
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and H P}Tn {z) a shorthand for the Hypergeometric func¬ 
tion i+i(— p/2, \m\ + 1, z). Relation (15) is obtained 
by inserting the expression of the CiB (2) into the in¬ 
tegral defining Ree(z), by changing variable into r = 
W(z)y/t/2, by taking the limit for z —► +oo, and by 
using the property lim 2 ^ +oc V2Ree(z)/W (z) = 9ee/9q- 


FIG. 2. Divergences 9 rms and See in function of p for the 
CiB with 3m(p) = 0, qi = 5Re(qo) and m = 2. When p is 
approaching to — |m|, the value S rms tends to infinity. 


To better understand the properties the CiBs, we now 
analyze their features for some particular values of the 
beam parameters and study in more detail their diver¬ 
gences and their experimental generation. Indeed, we 
would like to examine the CiBs with |£| = 1 correspond¬ 
ing to Qrn^i) = 0 or 3m(gi) —>• oo. In this case, a 
simpler expression for the CiB is obtained. Indeed, for 
|£| = 1 we have 2 = |p| 2 /[3?e(p) + \m\\ and the rms 

divergence and normalization factor simplify to 


9 r ms = 8q\ 1 + \m\ - 5Re(p^) + 


\p\ 2 Sfte(^) + 1 
2 3f ie(p) + \m\ 



T(|m| + 1 + 3?e(p)) 

|T(M + 1 + f)| 2 
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As it is evident from (17), when 3m (gi) = 0 and 3?e(p) = 
— \m\ the rms divergence cannot be defined, as already 
anticipated. For such beams, the divergence can be only 
evaluated by 9ee- To better illustrate the behavior of 
the divergences, we show in Figure 1 the values of 0 rms 
and 9ee in function of m for the CiB with qi = 3?e(go) ; 
such that £ = 1 and d rms = Ooy/l + m 2 /(p + |m|). We 
show three different CiBs classes, corresponding to the p 
parameter given by p = 0, p = — \m\ + 1 or p = —|m|. In 
the latter case (p = — |m|), as already discussed, only 9ee 
can be shown since 9 rms in not defined. The case p = 0 
corresponds to the LG mode LGo im . Similarly, in Figure 
2 we show the two divergences 9 rms and $ee i n function of 
p for the CiB with Sm(p) = 0, q\ = 3?e(go) and m = 2. 
It is worth noticing that increasing the value of p will 
reduce both divergences such that limp^+oo 9 = 9 0 . 

The role of the imaginary part of the p parameter can 
be appreciated when qi = —d\. In such case, at fixed 
5ie(p), the rms divergence is minimized for 3m (p) = 
(?fte(p) + |m|) and becomes a decreasing function 

of 3fte(p), namely 9 lms = 9 0 ^1 + ~J a+( lL d ^ ~ 
Then, when —q-[ = d\ ^ do, a complex value of p is 
required to minimize the beam divergence. 

We finally analyze a possible experimental generation 
of such subclass of CiBs. The cases qi = 3fJe(g 0 ) = — d 0 
and 3m(gi) —► oo respectively correspond to the HyGG 
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and HyGG-II modes introduced in [14, 15]. The general 
case qi = —di, that we call generalized Hypergeometric- 
Gaussian beam (gHyGG), is particularly interesting from 
the point of view of experimental generation. Indeed, the 
behavior of the field for z —> d± gives a clear hint for its 
possible experimental generation. For q\ = —d\ we have 

ikr 2 r p+\m\ e im<t> 

lim CiB^U^qx) oc e 2 (d 1 ~d 0 +iz 0 ) -- ; —. (18) 

z — >d,i ’ d i — do -h iZfl 

This feature is very useful because the r.h.s. of eq. (18) 
can be generated by applying the singular phase factor 
exp (im(f>) and a polynomial transmittance profile of the 
order p + |m| to Gaussian (TEM 0 o) beam at a distance 
di — do from its waist. In particular, the p = — |m| 
modes are simply generated by applying the phase factor 
exp {im<t>) to a Gaussian beam at a distance d\ — do from 
its waist plane. Such phase factor can be experimentally 
implemented by q-plates [ 20 ] or phase-plates [ 21 ]. 

In the present letter we studied the properties of the 
Circular-Beams, a general solution of the paraxial wave 
equation with OAM depending on three complex pa¬ 
rameters qo, <?i and p and an integer parameter to, re¬ 
lated to the content of OAM. The allowed values of the 


beam parameters are given in eq. ( 8 ) and the corre¬ 
sponding ones obtained by the symmetry (p, to, qo,qi) —> 
(—p — 2|to| — 2, to, < 71 , qo)- We derived the normalization 
of the CiBs and their expansion in terms of Laguerre- 
Gauss modes, a standard OAM basis (see eq. (4)). Such 
expansion allowed us to study the divergence of the CiBs 
in free-space propagation. We defined the divergence by 
two methods, the rnrs (see eq. (13)) and the encircled- 
energy (see eq. (15)) and studied their behavior for par¬ 
ticular values of the beam parameters. We finally sug¬ 
gested the experimental generation of a subclass of CiBs 
by looking at their behavior in a particular transverse 
plane, see (18). Our achievements, providing the esti¬ 
mation of the beam divergence for the CiBs (the most 
general beam with OAM known so far), may have po¬ 
tential application for OAM transportation in free space, 
such as multi-channel/multiplexing classical or quantum 
communication [ 6 , 10 , 22 ]. 

We thanks Paolo Villoresi of the University of Padova 
for suggestions and useful discussions. Our work was 
supported by the “Progetto di Ateneo PRAT 2013, OAM 
in free space: a new resource for QKD (CPDA138592)” 
of the University of Padova. 
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